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O ' Abstract 

<u : 

Q ■ A simple way to visualize event-by-event average pt fluctuations is by assuming that each collision 

has a different temperature parameter (inverse pt slope) and that the ensemble of events has a 
\ temperature distribution about the mean, (T), with standard deviation ot- PHENIX characterizes 

the non-random fluctuation of M PT , the event-by-event average pr, by F PT , the fractional difference 

O ' 

CN ■ of the standard deviation of the data from that of a random sample obtained with mixed events. 

This can be related to the temperature fluctuation: 

£ \ F PT = ojfc/o%£T - 1 * «n) - l)a 2 T /{Tf 



Combining this with I he Gavai. <! al.. 1} and Korus. el a£,Q dclinitjons of ihc spccilic Heal pvr 

particle, a simple relationship is obtained: 

3 _ Jn^^ 
^ ~(N tot )F PT ■ 

F PT is measured with a fraction (n) / {Ntot) of the total particles produced, a purely geometrical 
factor representing the fractional acceptance, ~ 1/33 in PHENIX. Gavai, et al. predict that c v /T 3 = 
15, which corresponds to F PT ~ 0.20% in PHENIX, which may be accessible by measurements of 
M PT in the range 0.2 < px < 0.6 GeV/c. In order to test the Gavai, et al. prediction that c v /T 3 is 
reduced in a QGP compared to the ideal gas value (15 compared to 21), precision measurements 
of F PT in the range 0.20% for 0.2 <pr< 0.6 GeV/c may be practical. 

PACS numbers: 24.60.-k, 25.75.-q, 25.75.Gz 
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I. INTRODUCTION 

R. Gavai, S. Gupta and S. Mukherjee predict in "quenched QCD" at 2T C and 3T C that 
the specific heat, cy/T 3 , differs significantly from the value for an ideal gas — 15 compared 
to 21 (see Fig. |TJ). Can this be measured? 
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FIG. 1: Gavai, et al, prediction for c v /T 3 
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II. EVENT-BY-EVENT AVERAGE p T FLUCTUATIONS AND SPECIFIC HEAT 
A. Single particle distributions 

The single particle transverse momentum (p?) distribution averaged over all particles in 
all events for a p-p experiment (inclusive) or in all events of a given centrality class for an 
A+A experiment (semi-inclusive) is usually written in the form: 

^ = ^(W H e^ or = JL( bpT r 2 e-^ . (1) 

Equation[T]represents a Gamma distribution, where (pr) = p/b, a PT / (pt) = 1/y/p- Typically 
6 = 6 (GeV/c)^ 1 and p = 2 for p-p collisions. As shown in Fig. El the p parameter depends 
on the particle type in central Au+Au collisions, with p < 2 for 7r ± , p ~ 2 for K ± and p > 2 
for (anti-) protons, but the asymptotic slope tends to be the same for all particles. The 
'inverse slope parameter', T = 1/6, is usually referred to as the 'Temperature parameter'. 
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FIG. 2: Identified particle semi-inclusive invariant pt spectra for Au+Au central collisions [3(]. At 
the lowest pr, the 7r + are the highest, followed by K + and p (left) and the same for the negatives 
(right). 

B. Event-by-Event Average 

For events with n detected charged particles with magnitudes of transverse momenta, p^ , 
the event-by-event average px, denoted M PT , is defined as: 
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By definition (M PT ) = (p?) = \i\ however, it takes hard work to make one's data follow this 
identity to high precision (<^C 1%). The standard deviation of M PT is defined the usual way: 



2 



PT 



M 



PT 



(M PT ) 



1 



li- 



no. 



Pt 



+ £ Z ((Pr* - (Pt))^ ~ (pr» 
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8=1 j=l,j^i 

If all the pr. on all events are random samples of the same distribution, then: 
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where <x 



PT 



(Pt) ~ (Pt) is the standard deviation of Eq. [TJ the inclusive px spectrum 



(averaged over all events). 

A nice illustration of what can be revealed by the event-by-event average that is not 

n 

shown by the inclusive average over all events was given by Korus, et al. |2||. Suppose that 
each collision has a different temperature parameter such that the ensemble of events has a 
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mean, (T), with standard deviation, ot = y (T 2 ) — (T) 2 , about the mean. It is easy to show 



that for this case: 



G kt-- a \ = ^--)^K • (5) 
fi 2 n fi 2 K n (T) 



C. Specific Heat 

ri 

As pointed out by Korus, et aZ., |2| if the parameter T would correspond to the actual 
temperature of the system, not just the inverse slope of the distribution, then a basic 
equation of thermodynamics would relate the temperature fluctuations of a system to its 
total heat capacity jj, 0, Q : 

Cy~ = Ixf ' (6) 

where Cy is an extensive quantity corresponding to the total number of particles in the 
system, (N tot ). Thus the specific heat per particle is cy = Cy j (N to t)- Gavai, et a/.,|l( refer 
to this same (dimensionless) quantity as c v /T 3 , resulting in the final equation: 

T 3 (N tot ) a 2 /(T) 2 1 ; 

where n represents the number of particles used in the calculation of M PT (Eq. EJ) from which 
o T j (T) is determined. 



III. MEASUREMENTS OF M PT 

The measured M PT distributions for two centrality classes in ^snn = 200 GeV Au+Au 
collisions in PHENIX j]| are shown in Figure El (data points) compared to a random baseline 
(histograms). Mixed-events are used to define the baseline for random fluctuations of M PT . 
This has the advantage of effectively removing any residual detector-dependent effects. The 
event-by-event average distributions are very sensitive to the number of tracks in the event 
(denoted n), so the mixed event sample is produced with the identical n distribution as the 
data. Additionally, no two tracks from the same data event are placed in the same mixed 
event in order to remove any intra-event correlations in pt- Finally, (M PT ) must exactly 
match the semi-inclusive (pt)- 

The non-Gaussian, Gamma distribution shape of the M PT distributions is evident. The 
difference between the data and the mixed-event random baseline distributions is barely 
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FIG. 3: M PT for 30-35% and 0-5% centrality classes|7|: data (points) mixed-events (histogram 
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visible to the naked eye. PHENIX quantifies the non-random fluctuation by the fractional 
difference of the standard deviations of M PT for the data and the mixed-event (random) 
sample: 



Pt 



®M Pt , data — 0"M Pt , mixed 
&M PT , mixed 



which is on the order of a few percent. The results are shown (Fig. EJ-left) as a 
function of centrality (represented by N part ) for charged particle tracks in the range 
0.2 GeV/c <p T < 2.0 GeV/c; and, for the 20-25% centrality class (N part = 181.6), over 
a varying p? range, 0.2 GeV/c < pt < p™ ax (Figure HI- right). The steep increase in F PT for 
the small increase in the number of tracks with increasing »^ ax > 1 GeV/c is consistent 



with correlations due to jet production as shown by the dotted lines |2|. However, other 
explanations have been proposed Note that the errors are entirely systematic, due to 
time-dependent detector variations. Comparatively, statistical errors are negligible. 



IV. HOW TO MEASURE c v /T 3 . 



For the small values of F PT observed, one can make use of the identity 
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FIG. 4: F PT vs centrality and p™ ax compared to simulations 
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to obtain the relation: 
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Combining Eq. fTTH with Eq. [?l we obtain the simple and elegant expression: 
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T 3 (Ntot) F PT 

Note that F PT is measured with a fraction (n) / (N to t) of the total particles produced, which 
is a purely geometrical factor representing the fractional acceptance of the measurement. 
For example, if all particles are produced in a range 5^fwhm (assuming a flat or trapezoidal 
dn/drj over this interval) and if including the neutrals gives a factor of 1.5 more total particles 
than charged particles; and if F PT is measured with charged particles in an acceptance Sr) c , 
S(J) c /2tt, which due to the pr cut only represents a fraction f c of the charged particles on that 
solid angle, then: 

1.5 x 2n x <5?7fwhm 



(n) f c x S(j) c x 5r/_ 

For RHIC at ^snn = 200 GeV, 5?7fwhm = ±3.5 and the PHENIX acceptance was 
5<p = 7T, 5r) = ±0.35, f c =0.9 for p T > 0.2 GeV/c, resulting in (N tot ) / (n) = 33. From 
Fig. El F PT is of order 2% but most of that is due to jets, so the effect due to temperature 
fluctuations is < F PT , say 1%, so we obtain c v /T 3 > 1/(33 * 0.01) = 3. This is to be 



compared to the Korus, et al. Q], result of c v /T'' (30 ± 100 



lOl from the NA49 data. Recall 



1] 



that Gavai, et ai, [l| predict a value of 15 for c v /T 3 , which would correspond to a value of 
F PT = 1% x 3/15 = 0.20% for the data in Fig. |U (see Fig. EJ. Perhaps this precision can 



c v /T 3 ~15 
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FIG. 5: Gavai et al., prediction [1| compared to PHENIX measurement. 



be achieved by concentrating on the region p™ ax < 0.6 GeV/c, where jets have least effect. 
Also, as the present error is totally systematic due to run-by-run variation, there is hope 
that a substantial reduction should be possible. 



V. CONCLUSIONS 

In central heavy ion collisions, the huge correlations in p-p collisions are washed out Q]. 
The remaining correlations are: Jets; Bose-Einstein correlations; Hydrodynamic Flow. These 
correlations seem to saturate the present fluctuation measurements. No other sources of 
non-random fluctuations have been observed. This puts a severe constraint on the critical 
fluctuations that were expected for a sharp phase transition but is consistent with the present 
expectation from lattice QCD that the transition is a smooth crossover. In order to see the 
temperature fluctuations predicted by c v /T 3 ~ 15 in lattice gauge calculations, present sen- 
sitivity needs to be improved by an order of magnitude by removing the known sources of 
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correlation and improving the measurement errors. An interesting check of whether tempera- 
ture fluctuations, rather than the correlations noted above, produce the observed non-random 
fluctuations is provided by Eq. [TUl for a pure a\j (T) 2 fluctuation, F PT for a given centrality 
should increase linearly with the number of tracks measured (e.g. by increasing the solid 
angle— PHENIX cf. STAR). 
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